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Abstract: soil physical properties (SPP) are considered to be important indices that reflect soil structure, 
hydrological conditions and soil quality. It is of substantial interest to study the spatial distribution of SPP 
owing to the high spatial variability caused by land consolidation under various land restoration modes in 
excavated farmland in the loess hilly area of China. In our study, three land restoration modes were 
selected including natural restoration land (NR), alfalfa land (AL) and maize land (ML). Soil texture 
composition, including the contents of clay, silt and sand, field capacity (FC), saturated conductivity (Ks) 
and bulk density (BD) were determined using a multifractal analysis. SPP were found to possess variable 
characteristics, although land consolidation destroyed the soil structure and decreased the spatial 
autocorrelation. Furthermore, SPP varied with land restoration and could be illustrated by the multifractal 
parameters of Dı, AD, Aa and Af in different modes of land restoration. Owing to multiple compaction 
from large machinery in the surface soil, soil particles were fine-grained and increased the spatial 
variability in soil texture composition under all the land restoration modes. Plough numbers and vegetative 
root characteristics had the most significant impacts on the improvement in SPP, which resulted in the 
best spatial distribution characteristics of SPP found in ML compared with those in AL and NR. In 
addition, compared with ML, Aa values of NR and AL were 4.9- and 3.0-fold that of FC, respectively, and 
Aa values of NR and AL were 2.3- and 1.5-fold higher than those of K,, respectively. These results 
indicate that SPP can be rapidly improved by increasing plough numbers and planting vegetation types 
after land consolidation. Thus, we conclude that ML is an optimal land restoration mode that results in 
favorable conditions to rapidly improve SPP. 
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1 Introduction 


Land consolidation is widespread throughout the world. Its goal is to provide a broad area 
of cultivation, improve the soil quality and production and ameliorate the ecological 
environment (Altes and Sang, 2011; Ying et al., 2020). There is an urgent need to consolidate 
the land in the 
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Loess Plateau in China, because this work is closely related to land development, agricultural 
production, soil desertification and flood control in the area of Yellow River (Zhu, 1995). 

A land consolidation project known as "Gullies Reclamation for Farmland" was implemented 
in the loess hilly area of China in 2013 (Ma et al., 2020). We designated the farmland formed after 
land consolidation as the excavated farmland, which is a large area of flat land formed by the 
large-scale mechanical excavation of slope soil landfill gully, after many times of compactions 
(Fig. 1). Land consolidation can easily cause spatial variability in soil physical properties (SPP) 
(Chen et al., 2015; Chen et al., 2019). The excavated slope soil has different soil textures owing to 
different ages of loess formation (Li et al., 2018), and the soil texture is unevenly distributed by 
random landfill. In addition, the filling of gullies with depths of several tens of meters may cause 
the groundwater level to deepen and uneven landfill and compaction lead to uneven soil bulk 
density (BD) and many void pores. These factors drastically affect the ability of water to infiltrate 
the soil and affect the field water holding capacity and even have extreme effects on the 
hydrological characteristics of the soil, including the surface water, groundwater and soil water 
(Tripathi et al., 2009; Marschalko et al., 2012; Ren et al., 2016). Moreover, an uneven SPP will 
cause spatial variation on water and nutrient storage capacity and even destroy agricultural 
ecosystems (Liang and Wei, 2020), which will seriously affect agronomic management of the 
fields. 

Studies have shown that land restoration can ameliorate SPP (Abed et al., 2020; Donovan and 
Monaghan, 2021; Wang et al., 2021), including a reduction in the rate of soil disintegration and 
BD (Wang et al., 2017), an increase in the soil moisture content (Fu et al., 2020), and enhanced 
conductivity of the saturated water and the field water holding capacity (Li et al., 2017; Wang et 
al., 2017; Mehdi et al., 2019; Lozano-Baez et al., 2019). In addition, Perring et al. (2012) showed 
that the restoration of agricultural land can benefit the ecosystem in multiple manners. To this end, 
three land restoration modes that included natural restoration (NR), alfalfa (Medicago sativa Linn) 
land (AL) and maize (Zea mays L.) land (ML) were implemented in excavated farmland (Fig. 2). 
Based on natural restoration of abandoned farmland in the Loess Plateau, SPP directly related to 
the vegetation recovery stages and the time of abandonment. Rate of improvement of the physical 
property of 0-20 cm soil layer is significantly higher than that of 20-40 cm soil layer (Li and 
Shao, 2006). Dong et al. (2016) reported that there was no considerable variation in the soil 
properties between AL and common farmland 10 years ago. However, after 10 years, the soil 
water holding capacity and water stability were significantly higher than that of common 
farmland, indicating that AL improved soil environment more favorably. It is apparent that the 
types of vegetation and tillage years have various effects on SPP. Therefore, we assumed that 
there are some variations in the spatial distribution of SPP under different land restoration modes 
in the excavated farmland. 

To obtain the spatial variability of SPP, many researchers interpret spatiotemporal distributions 
at different scales and use appropriate research methods for data processing (Paterson et al., 2018; 
Li et al., 2019; Mojtaba et al., 2019). Classical statistics describe the level of variations of SPP 
using a coefficient of variation, but this method seriously ignores the relative position between 
variables. Geostatistics can effectively explain the spatiotemporal autocorrelation of variables 
using a semi-variant function and Kriging, but it is difficult to characterize the spatial variation 
(Premo, 2004). Traditional methods have an advantage since they use simple and independent 
analyses, but practically there are several complicated parameters that are not fully independent 
(Jing et al., 2019). Fractal theory based on independent correlation of variables is widely used in 
soil sciences, meteorology and information sciences to solve various complicated problems at the 
micro and macro levels (Morató et al., 2017; Paterson et al., 2018; Wang et al., 2018; Gao et al., 
2021; Stanić et al., 2021; Wu et al., 2021). Qi et al. (2018) and Xia et al. (2020) reported that the 
spatial distribution difference of soil texture under different vegetation types can be effectively 
characterized using multifractal theory. Jing et al. (2019) showed that BD, saturated water 
conductivity and water holding capacity, all have good multifractal characteristics in farmland. 
Therefore, multifractal theory is an appropriate choice to evaluate the spatial variability of SPP. 

The spatial distribution of SPP not only determines the soil water, fertilizer, gas and heat, but 
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also affects the availability and supply of plant nutrients in the soil (Drewry, 2006). Therefore, it 
is considered to be an important index to reflect soil structure and hydrological conditions, and to 
evaluate soil quality (Doran et al., 1996; Boix-Fayos et al., 2001). However, there are still gaps in 
the study of spatial distribution of SPP under different land restoration modes in the excavated 
farmland due to sustainable land consolidation. Consequently, the purposes of this study are (1) to 
characterize the spatial variability of SPP using multifractal analysis in the excavated farmland; 
and (2) to discuss the mechanism of the influence of different land restoration modes on SPP. 


2 Materials and methods 


2.1 Study area 


The study area is located in an alluvial farmland of loess hilly region in China (40°14'11"N 
—42°27'42'"'N, 75°33'16'"E—80°59'7"E; Fig. 1). The region has a semi-arid continental monsoon 
climate with a mean annual precipitation of 505 mm and mean annual potential evaporation of 
1463 mm (Ke et al., 2021). 
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Fig. 1 (a), location of Ansai District in the hilly loess region of China; (b), schematic illustration of the 
"Gullies Reclamation for Farmland" project; (c), excavated farmland. 


For this study, data were collected from three land restoration modes that included NR, AL and 
ML in the excavated farmland (Fig. 2). The topography of NR, AL and ML was similar. Moreover, 
one NR, one AL and one ML plot were selected and all the plots were restored after 5 years. The 
plot sizes were kept the same, i.e., 60 mx60 m. 


Fig. 2 Overview of study area. (a), natural restoration land; (b), alfalfa land; (c), maize land. 
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2.2 Agronomic management 


NR was undisturbed, and its vegetation type is mostly annual herbaceous plants. AL was 
five-year-old alfalfa. The forage was harvested at three times a year and ploughed once a year 
before returning to green. ML was planted to a season of spring maize, which was ploughed and 
weeded twice during the growth season, and total 600 kg P2Os/hm? (diammonium phosphate) and 
100 kg N/hm? (urea) were applied as base fertilizer. 


2.3 Sample collection and measurement 


A sampling network along with squared grids (15 mx15 m) was constructed in each 0.36 hm? plot. 
NR, AL and ML plots consisted of 48 sampling points. Each sampling point was located at the 
individual grid center. Soil profiles were carried out in July 2018. Trimble GeoXT GPS 
equipment (Sunnyvale, CA, USA) with a positioning precision <0.5 m was used to locate the 
sampling center point coordinates. Soil texture (clay, silt and sand content), BD, saturated 
hydraulic conductivity (Ks) and field capacity (FC) were measured in 0—20 cm and 20—40 cm soil 
layers. At each sampling point, undisturbed soil samples were analyzed for their BD, FC and Ks 
using a ring knife (100 cm). To examine the soil particle size distribution (PSD), we first 
removed the roots and impurities in air-dried soil samples and passed through a 2-mm screen, 
before adding a solution of 10% H2O2 to degrade any organic matter. Then 0.3 g of soil was 
employed for laser diffraction analysis (Mastersizer 2000, Malvern Company, UK). BD and FC 
were evaluated using cutting ring knife method (Federer, 1983; Prévost, 2004). Ks was evaluated 
by a constant elevation head permeameter (TST-55, Beijing Aerospace Huayu Test Instrument Co., 
Ltd., Beijing, China) (Li and Shao, 2006). 


2.4 Calculation of multifractal parameters 


In this study, spatial variability of SPP was evaluated by multifractal analysis. A grid square with 
size €, enclosed by a part of space that has SPP. These parameters were selected to identify 
singularity exponent a(q), generalized dimension D(q) and singularity spectrum f(a) (Jing et al., 
2019). In each plot, there were N(e)=2* (k=0, 1, 2, ...) cells that were chosen (Caniego et al., 2005; 
Morató et al., 2017). Here in this paper, a scale of 60 mx60 m was used for each sampling plot 
while a scale of 15 mx15 m was used as a sampling point size. In this work, we considered four 
grid sizes (15, 20, 30 and 60 m) with a total number of grids, i.e., 16, 9, 4 and 1 in each plot 
accordingly. 

A measure was defined with respect to the underlying steps to conduct a multifractal analysis. 
First, P(¢) named the probability mass function was calculated as follows: 


Z, 
a en (1) 
E 
where Z; is the value of the measure in a given size e, and N(e) is the number of grids. 
Secondly, generalized fractal dimension (D,) is defined as follows: 
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where q is the integer in (+œ, —), i.e., the probability density weight index (Li et al., 2011). In 
this study, generalized fractal dimensions D, of reconstructed SPP were calculated for —10<q<10 with 
increment of 1. D; is the information dimension. A relatively lower value of Dı showed that the value 
of SPP was distributed over a relatively larger domain, owning relatively higher spatial variability, 
while a relatively higher value of Dı showed that the value of SPP was focused on a relatively smaller 
domain, owning relatively lower value of spatial variability (Evertsz and Mandelbrot, 1992; Jing et al., 
2019). AD, D-1o—D1o, can evaluate the intensity of spatial variability of SPP in a local distribution. AD 
value is directly correlated with variability (Zhou et al., 2010). 
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Thirdly, measure's multifractal spectrum, fq), can express the results. The following Legendre 
transformation defines fiq): 
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where u:(q, €) is the probability of q at the subinterval i; and a(q) is the singularity exponent. a(q) 
and f(q) were calculated by using Equations 4 and 5. Multifractal spectrum width A@ (Qmax—Qmin) 
can reflect the overall variability in spatial distribution. Af, f(@min)—/(Qmax), reflects the 
multifractal spectrum shape feature and it characterized the degree of symmetry for SPP. Af<0 
shows that a subset is dominant with a small probability and f(g) has a form similar to that of a 
hook to the right. In contrast, Af>0 denotes that a subset is dominant with a large probability and 
J(q) has a form similar to a hook to the left (Morató et al., 2017). 

2.5 Statistical analysis 

The values calculated for SPP, i.e., Dı, AD, Aa and Af, were evaluated using statistical methods. 
SPSS v23.0 (IBM, Inc., Armonk, NY, USA) was used to determine the mean, coefficient of 
variation (CV) and standard deviation (SD). Least significant difference (LSD) tests were 
significant at P<0.05 level. CV was ranked into three levels that were smaller (<10%), middle 
(10%-100%) and greater (>100%). Graphs were prepared using Origin Pro v8.0 (OriginLab, 
Northampton, MA, USA). 


3 Results 


3.1 Descriptive statistics of SPP 

Table 1 showed the mean values, CV and significant difference of SPP in 0-20 cm and 20-40 cm 
soil layers under three modes of land restoration. The order of clay contents was ML>AL>NR in 
the whole layer of soil (P<0.05), and sand contents were ML<AL<NR (P<0.05). CVs of clay, silt 
content and sand were all lower than <0.10 in the whole soil layer, indicating a low variability. 


Table 1 Soil bulk density (BD), field capacity (FC), saturated hydraulic conductivity (Ks) and soil texture (clay, 
silt and sand contents) in three plots 


Soil property iby ase in ML CV AL CV NR CV 
Clay content 0-20 9.85+0.34^ 0.03 9,6240.364 0.04 9.490.364 0.04 
(%) 20-40 9,990.32" 0.03 9.7240.314 0.03 9.61+0.36% 0.04 

Silt content 0-20 39.20+1.72"" 0.04 39,33+1.74" 0.04 38.61+1.82^ 0.05 
(%) 20-40 40.04+1.71% 0.04 39.36+1.7™ 0.04 38.97+1.70* 0.04 
Saideontent 020 50.95+1.61*" 0.03 51.06+1.57% 0.03 51.89+1.69^ 0.03 
(%) 20-40 49.97+1.69^ 0.03 50.92+1.68^? 0.03 51.42+1.68^ 0.03 
0-20 1.06+0.088° 0.07 1.31+0.078° 0.05 1.52+0.07^ 0.05 

BD (gom 49-40 1.2540.094° 0.07 1.4940.114° 0.07 1.55£0.12 0.08 
0-20 24.26+0.81^ 0.03 21.44+1.29^ 0.06 19.48+1.58^° 0.08 

FC %) 20-40 23.21+0.7288 0.03 20.08+1.108° 0.05 18.91+1.494¢ 0.08 

l 0-20 0.6240.13" 0.20 0.47+0.24%° 0.50 0.350.254 0.71 

Ka (emmin) aa 40 0.39+0.17®° 0.44 0.20+0.12®® 0.58 0.16+0.16®” 0.99 


Note: ML, maize land; AL, alfalfa land; NR, natural recovery land; CV, coefficient of variation. Different uppercase letters within the 
same plot indicate significant differences between depths at P<0.05 level. Different lowercase letters within the same depth indicate 
significant differences among plots at P<0.05 level. Mean+SD. 
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In the whole soil layer, BD of NR were 32.9% and 9.6% and significantly higher than those of 
ML and AL, respectively (P<0.05). FC of ML were 14.3% and 23.7% and significantly higher 
than those of AL and NR, respectively (P<0.05). Ks of ML were 50.7% and 98.0% and 
significantly higher than those of AL and NR, respectively (P<0.05). BD and FC of all the land 
restoration modes had a low degree of variability over the whole soil layer. Ks of ML had a 
middle variability in 0-20 cm soil layer, while the others were highly variable. 


3.2 Multifractal analysis of SPP 


3.2.1 Generalized dimension of spatial distribution 
Generalized dimension spectra of SPP in 0-20 cm and 20—40 cm soil layers under ML, AL and 
NR are shown in Figures 3 and 4. D, values inversely decreased with q values as shown in the 
spectra "S" shape curve, which indicate that the variables have multifractal characteristics. In this 
study, SPP having multifractal characteristics implied that multifractal evaluation could be 
applied. 

Table 2 shows the multifractal parameters of SPP. As for the soil texture, D; showed the 
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Fig. 3 Generalized fractal dimension (D4) of soil texture (clay (a and b), silt (c and d) and sand contents (e and 
f)) in different soil layers in three plots. g, probability density weight index; ML, maize land; AL, alfalfa land; NR, 
natural recovery land. 
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Fig. 4 Generalized fractal dimension (D4) of soil bulk density (BD; a and b), field capacity (FC; c and d) and 
saturated hydraulic conductivity (Ks; e and f) in different soil layers in three plots. g, probability density weight 
index; ML, maize land; AL, alfalfa land; NR, natural recovery land. 


tendency of ML>AL=NR over the whole soil layer and whatever clay, silt and sand contents are. 
D; values were higher in 20—40 cm layer than in 0-20 cm layer under the same mode, and vice 
versus for AD. For BD, Dı showed the tendency of ML<AL<NR over the whole layer of soil, 
while AD was ML>AL>NR. For FC and Ks, Dı and AD showed the similar tendency with soil 
texture in different land restoration modes. Dı value was higher in 20-40 cm layer than in 0-20 
cm in FC, which was opposite for Ks. AD was lower in 20—40 cm layer than in 0-20 cm in FC, 
which was opposite for Ks. 
3.2.2 Singularity spectra of SPP 
The singularity spectra of SPP was shown in Figures 5 and 6. Aa values of clay, silt and sand were 
higher in 0-20 cm layer than in 20—40 cm soil layer under all the land restoration modes. For BD, 
Aa showed the tendency of ML>AL>NR in the whole soil layer, which was opposite for FC and 
Ks. Compared with ML, Aa values of NR and AL were 4.9- and 3.0-fold that of FC, respectively, 
while the comparable values were 1.5- and 2.3-fold that of Ks. 

Af values of clay and sand were greater than 0, whose curve was a left hook under all the land 
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Table 2 Multifractal parameters for soil bulk density (BD), field capacity (FC), saturated hydraulic conductivity 
(Ks) and soil texture (clay, silt and sand contents) in three plots 


Soil Soil depth ML AL NR 
property (cm) D AD ha Af D AD da Af D AD da Af 
Clay content 0-20 9992 0.0169 0.0335 0.0081 1.9990 0.0195 0.0388 0.0091 1.9990 0.0198 0.0395 0.0229 
20-40 1.9994 0.0125 0.0253 0.0225 1.9994 0.0124 0.0249 0.0222 1.9992 0.0168 0.0335 0.0229 
Mean 1.9993 0.0147 0.0294 0.0153 1.9992 0.016 0.0319 0.0157 1.9991 0.0183 0.0365 0.0229 
Silt content 0-20 9988 0.0238 0.0469 -0.0296 1.9988 0.0238 0.0470 -0.0272 1.9986 0.0279 0.0551 -0.0355 
20-40 1.9993 0.0146 0.0288 -0.0245 1.9992 0.0158 0.0312 -0.0259 1.9992 0.0163 0.0320 -0.0269 
Mean 1.9991 0.0192 0.0379 -0.0271 1.9990 0.0198 0.0391 -0.0266 1.9989 0.0221 0.0436 -0.0312 
Sand content 0-20 9994 0.0122 0.0245 0.0111 1.9994 0.0117 0.0235 0.0103 1.9993 0.0133 0.0267 0.0129 
20-40 1.9995 0.0101 0.0204 0.0187 1.9995 0.0101 0.0203 0.0163 1.9995 0.0095 0.0190 0.0154 
Mean 1.9995 0.0112 0.0225 0.0149 1.9995 0.0109 0.0219 0.0133 1.9994 0.0114 0.0229 0.0142 
BD 0-20 9965 0.0688 0.1356 0.1243 1.9987 0.0246 0.0476 -0.0440 1.9990 0.0199 0.0390 -0.0096 
20-40 1.9969 0.0596 0.1150 0.1126 1.9972 0.0546 0.1066 -0.0138 1.9974 0.0518 0.1036 0.1345 
Mean 1.9967 0.0642 0.1253 0.1185 1.9980 0.0396 0.0771 0.0289 1.9982 0.0359 0.0713 0.0625 
FC 0-20 9994 0.0126 0.0249 -0.0041 1.9976 0.0462 0.0890 0.0256 1.9967 0.0614 0.1134 -0.0257 
20-40 1.9995 0.0103 0.0204 0.0018 1.9987 0.0253 0.0488 -0.0119 1.9974 0.0539 0.1096 0.0338 
Mean 1.9995 0.0115 0.0227 0.0012 1.9982 0.0358 0.0689 0.0069 1.9971 0.0577 0.1115 0.0041 
K, 0-20 9745 0.3820 0.6086 -0.2603 1.8837 1.0203 1.2591 0.7595 1.7287 1.6684 1.9590 1.3434 
20-40 1.9047 1.0157 1.3356 0.5495 1.7815 1.5097 1.7385 -0.7527 1.6603 2.1597 2.4938 -0.7405 
Mean 1.9396 0.6989 0.9721 0.1446 1.8326 1.2650 1.4988 0.0034 1.6945 1.9141 2.2264 0.3015 


Note: Dı, information dimension; AD, D_,o—Djo, spatial variability in a local distribution; A@, Q@max—Q@min, Spatial variability in an overall 
distribution; Af, A admin) -f max), multifractal spectrum shape feature; ML, maize land; AL, alfalfa land; NR, natural recovery land. 
restoration modes, whereas silt was a right hook (Fig. 5). For BD, the spectral shapes of AL and 
NR revealed a right hook in contrast to ML in 0-20 cm soil layer, and the spectral shapes of ML 
and NR revealed a right hook in contrast to AL in 20—40 cm soil layer (Fig. 6). For FC, spectral 
shape of ML and NR revealed a right hook and AL revealed a left hook in 0-20 cm soil layer in 
contrasted to 20—40 cm soil layer. For Ks, the spectral shapes of ML, AL and NR revealed a left 
hook in the whole soil layer. 


4 Discussion 


4.1 Effects of land restoration modes on SPP 


The homogeneity of SPP seriously affects soil moisture and nutrient distribution, which is an 
important consideration in farmland management (He et al., 2019). A poor soil structure and a high 
spatial variability of soil properties after land consolidation have been reported by Chen et al. (2015, 
2019). Therefore, many land restoration modes were used to ameliorate soil properties to facilitate 
crop growth and increase yield (Lozano-Baez et al., 2019; Dou et al., 2020; Sipek et al., 2020). 

In our study, the composition of soil texture changed slightly, and BD, FC and K, were 
significantly different under the three modes after 5 years of land restoration (Table 1). We found 
that the spatial variabilities of soil texture were higher in 0—20 cm soil layer than in 20—40 cm soil 
layer under all the land restoration modes. Moreover, the contents of clay and silt were higher in 
0-20 cm soil layer than in 20—40 cm soil layer, in contrast to sand content (Table 2). These results 
indicated that the content of small size of soil particles increased, while the distribution 
uniformity decreased in 0-20 cm soil layer. Min et al. (2017) investigated PSD in reclaimed soil 
produced by mechanical compaction and found that soil particle sizes became finer with the 
number of compactions. In our study, surface soil was compacted and leveled several times 
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Fig. 5 Multifractal spectrum functions for soil texture (clay (a and b), silt (c and d) and sand contents (e and f)) 
in different soil layers in three plots. f(a), singularity spectrum; a(q), singularity exponent; ML, maize land; AL, 
alfalfa land; NR, natural recovery land. 


during the formation of excavated farmland, suggesting that a large part of soil particles changed 
into fine particles and increased the spatial variability of soil texture in 0-20 cm soil layer. 
Moreover, there were significant differences in the compositions of soil texture under different 
land restoration modes (Table 1). A difference in the vegetation cover can be a possible primary 
reason. Qi et al. (2018) have shown that soils in oak forestland, shrub grass sloping land, terraced 
farmland and sloping farmland were classified as silty loam, sandy loam, sandy loam and loamy 
sand, respectively, in the Funiu mountainous region, China. Liu et al. (2009) have also shown 
considerable differences in PSD among seven plant communities in the Yimeng Mountain, China. 
Another reason of the changes in soil texture could be the frequency of mechanical disturbance. 
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Fig. 6 Multifractal spectrum functions for soil bulk density (BD; a and b), field capacity (FC; c and d) and 
saturated hydraulic conductivity (Ks; e and f) in different soil layers in three plots. f(a), singularity spectrum; 
a(q), singularity exponent; ML, maize land; AL, alfalfa land; NR, natural recovery land. 


ML was ploughed and weeded twice during growth season. Forage of AL was mechanically 
harvested three times a year and once a year before returning green. NR was not disturbed by 
mechinary. Summary, the frequency of mechanical disturbance was the main factor that 
influenced the changes in soil particles, which was consistent with the result of Min et al. (2017). 

ML showed the most significant impacts on BD, FC and Ks, followed by AL and NR (Tables 1 
and 2). There are two explanations that consist of both plough numbers and vegetation root 
characteristics that help explain this result. (1) Owing to a lack of human disturbance, NR 
basically maintained the condition in which the soil was relatively uniformly compacted after the 
farmland creation. However, since AL was ploughed once a year before returning to green, and 
ML was planted to a season of spring maize, and ploughed and weeded twice during growth 
season, which indicated that, with the increase in plough times, BD decreased, while FC and Ks 
increased. Thus, plough times were the primary reason that caused ML to be more significantly 
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impacted on BD, FC and K, than those of AL. (2) Most of vegetation types were annual herbs in 
NR. Guo et al. (2018) reported that the root density of abandoned land was 1.54 kg/m? and their 
roots were concentrated in 0-10 soil layer after 8 years of natural recovery. Huang et al. (2019) 
showed that root densities of alfalfa and maize were 13.1 and 6.0 kg/m’, respectively, and 80% of 
the roots were concentrated in 0-30 cm soil layer. Soil porosity increased in the area with dense 
root distribution. Therefore, compared with NR, AL and ML had more significant effects on the 
values and spatial variability of BD, FC and Ks. Moreover, plough numbers were considered 
dominant because ML had a better SPP than that of AL. These results implied that SPP can be 
rapidly improved by increasing plough numbers and planting vegetation types with a high root 
density after land consolidation. In our study, we found that values of BD, FC and K, were close 
to those of sloping farmland and terraces in this area after 5 years of land restoration (Wang et al., 
2008), thereby indicating that this mode could allow crops to become as productive as normal. In 
summary, we conclude that ML is an optimal land restoration mode that results in favorable 
conditions for the rapid improvement of SPP. 


4.2 Application of multifractal analysis in SPP 


Spatial variability analysis of SPP plays a significant hypothetical role in understanding the 
excavation and management of farmland. In this study, we observed that physical characteristics 
of soil have multifractal characteristics using multifractal analyses and parameters (1.e., Di, AD, 
Aa and Af), which explained the concentration, local variability, overall variability and 
distribution symmetry of SPP in detail (Figs. 3—6). Using a traditional statistical analysis, we 
found that CVs of soil texture compositions, BD and FC were in the range of 0.03—0.08. CVs of 
both soil layers (0-20 and 20-40 cm) were similar, but that does not mean that their spatial 
distribution was similar (Table 1). Teng et al. (2017) showed that the spatial autocorrelation of 
soil organic carbon can be effectively expressed through geostatistical analysis, while description 
of spatial variation was uneven. Studies have shown that a fractal dimension can only reflect the 
spatial variation but cannot describe the characteristics of PSD more carefully in various types of 
vegetation using a single fractal analysis (Liu et al., 2009). Moreover, Liao et al. (2017) also 
found limitations in studying spatial and temporal patterns of soil water content. These results 
indicated that a multifractal analysis that could describe spatial variations in more detail and in 
many dimensions and had a substantial advantage in studying the spatial variation of SPP. 

Multifractal analysis was based on autocorrelation that was an inherent characteristic of soil 
properties (Morató et al., 2017). Qi et al. (2018) and Wang et al. (2018a) showed that multifractal 
analysis can effectively characterize the traits of distribution of soil particle size. Zhang et al. 
(2019) analyzed soil moisture and soil particle size on the slope of Loess Plateau using a 
multifractal analysis and found that there was a close relationship between them. We found that 
most of studies used multifractal analyses to study natural soil, and few of them studied soil 
properties after land consolidation. Land consolidation destroys soil structure and substantially 
reduces spatial autocorrelation of SPP. In addition, Pachepsky and Kravchenko (2004) delineated 
that this type of soil has no multifractal characteristics in some properties. In our study, ML and 
AL that had been disturbed by large machinery exhibited good multifractal characteristics, and 
this was also the case for NR. Jing et al. (2019) showed that SPP in subsidence areas in coal 
mines had multifractal characteristics, and changes in the effects of soil characteristics between 
land settlement and land restoration were effectively described by multifractal parameters. 
Therefore, we concluded that a multifractal analysis can be used to study SPP in more 
complicated situations. 


5 Conclusions 


In this study, spatial variability of SPP (sand, clay, silt, BD, FC and Ks) were analyzed, and 
mechanism of effects of diverse land restoration modes on SPP were illustrated using a 
multifractal analysis for three land restoration modes (ML, AL and NR) after 5 years in the 
excavated farmland. The main conclusions are as follows: under all the land restoration modes, 
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soil particles had fine grains owing to multiple compactions of large machinery in 0—20 cm soil 
layer. Plough numbers and vegetation root characteristics showed the most significant impacts on 
the improvement of SPP. These results indicated that ML exhibited the optimal spatial distribution 
characteristics of SPP compared with those of AL and NR. In addition, values of BD, FC and Ks 
for ML were close to those of common farmland in this area after 5 years of land restoration. 
These results implied that ML was an optimal land restoration mode that resulted in favorable 
conditions for the rapid improvement of SPP. 

Land consolidation destroyed soil structure and substantially reduced spatial autocorrelation of 
SPP. However, ML, AL and NR still had multifractal characteristics. Moreover, multifractal 
parameters (i.e., Di, AD, Aa and Af) explained the concentration, local variability, overall 
variability and distribution symmetry of SPP in detail. Therefore, we concluded that multifractal 
analysis can be used to study SPP in more complicated situations. 
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